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Abstract 

Let V he a system of n linear nonhomogeneous ordinary differential polyno- 
mials in a set ?7 of n — 1 differential indeterminates. Differential resultant 
formulas are presented to eliminate the differential indeterminates in U from 
V. These formulas are determinants of coefficient matrices of appropriate 
sets of derivatives of the differential polynomials in V, or in a linear pertur- 
bation Vs of V. In particular, the formula 9FRes('P) is the determinant of a 
matrix Ai{V) having no zero columns if the system V is "super essential". 
As an application, if the system ^ is sparse generic, such formulas can be 
used to compute the differential resultant 9Res(*^) introduced by Li, Gao 
and Yuan in [20]. 
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differential resultant, linear perturbation 
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1. Introduction 

Elimination theory has proven to be a relevant tool in (differential) al- 
gebraic geometry (see [8], [9] and [1]). Elimination techniques have been de- 
veloped using Grobner bases, characteristic sets and (differential) resultants. 
The algebraic resultant has been broadly studied, regarding theory and com- 
putation, some significant references are [14], [5], [28] and [11]. Meanwhile, 
its counterpart the differential resultant is at an initial state of development, 
a survey on this development can be found in the introductions of [15] and 
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[25]. Until very recently, the existing definitions of differential resultants for 
differential polynomials depended on the computation method [6]. In the 
recent paper [15], a rigorous definition of the differential resultant 9Res(*P), 
of a set *^ of n nonhomogeneous generic ordinary differential polynomials in 
n — 1 differential variables, has been presented: If the elimination ideal, of 
the differential ideal generated by is — 1 dimensional then it equals the 
saturation ideal of a differential polynomial 9Res(*P), the differential resul- 
tant of As in the algebraic case, the object that is naturally necessary for 
applications is the sparse differential resultant, and this was defined in [20], 
for a set of nonhomogeneous generic sparse ordinary differential polynomials. 

The computation and applicability of sparse algebraic resultants attained 
great benefits from having close formulas for their representation [11], [29], 
[12]. These formulas provide bounds for the degree of the elimination output 
and ways of exploiting sparseness of the input polynomials on predicting the 
support of the elimination output. Namely, obtaining the Newton polytope of 
the resultant [28], whose support is a superset of the support of the resultant, 
reduces elimination to an interpolation problem in (numerical) linear algebra, 
[29], [10], [12]. 

Sparse differential resultants can be computed with characteristic set 
methods for differential polynomials via symbolic computation algorithms 
[4], [18], [15], [23]. The algorithms in [18] and [4] have been implemented 
in the Maple package diffalg, [3] and in the BLAD libraries [2] respectively. 
These methods do not have an elementary complexity bound [17] but, a sin- 
gle exponential algorithm based on order and degree bounds of the sparse 
differential resultant has been recently proposed in [21]. It would be useful 
to represent the sparse differential resultant as the quotient of two deter- 
minants, as done for the algebraic case in [11]. As noted in [20] and [21], 
having similar formulas in the differential case would improve the existing 
bounds for degree and order of the sparse differential resultant and therefore 
the existing algorithms for its computation. Matrix formulas would also con- 
tribute to the development of methods to predict the support of the sparse 
differential resultant, achieving similar benefits to the ones obtained in the 
algebraic case. A matrix representation of the sparse differential resultant 
is important because it is the basis for efficient computation algorithms and 
their study promises to have a grate contribution to the development and 
applicability of differential elimination techniques. 

In the differential case, these so called Macaulay style formulas do not 
exist, even in the simplest situation. The matrices used in the algebraic case 
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to define the Macaulay style formulas [11], are coefficient matrices of sets of 
polynomials obtained by multiplying the original ones by appropriate sets of 
monomials, [5]. In the differential case, in addition, derivatives of the original 
polynomials should be considered. The differential resultant formula defined 
by Carra-Ferro in [6], is the algebraic resultant of Macaulay [22], of a set of 
derivatives of the ordinary differential polynomials in V. Already in the linear 
sparse generic case, these formulas vanish often, giving no information about 
the differential resultant 9Res(^). A determinantal formula for 2 generic 
differential polynomials of arbitrary degree and order 1 has been recently 
presented in [30]. 

In this paper, determinantal formulas are provided for systems of n linear 
nonhomogeneous ( non necessarily generic) differential polynomials P in n— 1 
differential indeterminates. The linear case can be seen as a previous study to 
get ready to approach the nonlinear case. One can consider only the problem 
of taking the appropriate set of derivatives of the elements in V and forget 
about the multiplication by sets of monomials for the moment. 

Given n differential polynomials, differential elimination is guaranteed of 
at most n — 1 differential variables (see Section 5) but, if there were more, 
we may decide which ones to consider as part of the coefficients. Take for 
instance the Lotka-Volterra equations 

J x' = ax — I3xy, 

1 y' = iy- pxy, 

they can be looked at as a system given by two linear differential polynomials 
in the differential indeterminate x, with a, /3, 7 and p algebraic constants, 

fi{x) = {f3y — a)x + x' = aix + a2x', 
/2(x) =y' --fy + pyx = bo + biX. 

Elimination of the x differential variable can be achieved by the determinant 
of the coefficient matrix of fi{x), f2{x) and f2{x), 

diiyT - yy" + ^-yy' - acy^ - hy'^y' + hcy^). 

In [27], the linear complete differential resultant dCKes{V) of a set of 
linear differential polynomials V was defined, as an improvement, in the 
linear case, of the differential resultant formula given by Carra-Ferro. Still, 
dGKes{V) is the determinant of a matrix having zero columns in many cases. 
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An implementation of the differential resultant formulas defined by Carra- 
Ferro and the linear complete differential resultant defined in [26] is available 
at [24]. 

The linear differential polynomials in V can be described via differential 
operators. We use appropriate bounds of the supports of those differential 
operators to decide on a convenient set ps(P) of derivatives of V, such that its 
coefficient matrix A4.{V) is squared and, if V is super essential (as defined in 
Section 3), it has no zero columns. Obviously, det{M.{V)) could still be zero. 
In such case, we can always provide a linear perturbation Ve of V such that 
det^M-iVe)) 7^ 0, as an adaptation of the perturbation methods described 
in [26] (for linear complete differential resultants) to the new formulas pre- 
sented in this paper. In the sparse generic case, we can guarantee that the 
linear sparse differential resultant 9Res(*^) can always be computed via the 
determinant of the coefficient matrix Aii^l) of a set ps(*^*), of derivatives 
of the elements in the perturbation of a super essential subsystem of 

Given a system of linear nonhomogeneous ordinary differential polyno- 
mials P, in Section 2, we describe appropriate sets bounding the supports 
of the differential operators describing the polynomials in V. Differential 
resultant formulas for V are given in Section 3. In particular, the formula 
9FRes('P) is defined, for the so called super essential (irredundant) systems, 
as the determinant of a matrix Ai(V) with no zero columns. In Section 4, it 
is shown that every system V contains a super essential subsystem V*, which 
is unique if V is differentially essential. Results on differential elimination 
for systems V of linear differential polynomial parametric equations (linear 
DPPEs) are given in Section 5, including a perturbation Ve such that if V 
is super essential then 9FRes('Pe) ^ 0. The methods in Section 5 are used 
in Section 6 to compute the differential resultant 9Res(*P) of a linear non- 
homogeneous generic sparse system ^ of ordinary differential polynomials. 
As explained in Section 6, the differential resultant 9Res(*P) exists only for 
differentially essential systems. 

2. Preliminary notions 

Let D be an ordinary differential domain with derivation d. Let us con- 
sider the set U = {ui, . . . , Wn-i} of differential indeterminates over D. By Nq 
we mean the natural numbers including 0. For k G No, we denote by Uj^k the 
kth derivative of Uj and for uj^ we simply write uj. We denote by {f/} the 
set of derivatives of the elements of U, {U} = {d^u \ u E U, k E No}, and by 
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D{[/} the ring of differential polynomials in the differential indeterminates 
U, which is a differential ring with derivation d, 



Given a subset U C {U}, we denote by D[W] the ring of polynomials in the 
indeterminates U. Given / G D{f/} and ?/ G f/, we denote by ord(/, y) the 
order of / in the variable y. If / does not have a term in y then we define 
ord(/, ?/) = —1. The order of / equals max{ord(/, ?/) | y G U}. 

Let /C be a differential field of characteristic zero with derivation d (e.g. 
/C = Q(t), d = d/dt) and C = {ci, . . . , c„} a set of differential indeterminates 
over /C. The differential ring /C{C} is an example of differential domain. By 
/C(C) we denote the differential field extension of /C by C, the quotient field 
of /C{C}. The following rankings will be used throughout the paper (see [19], 
page 75): 

• The order Ui < ■ ■ ■ < u^-i induces an orderly ranking on U (i.e. an 
order on {U}) as follows: Uij < u^^i <^ (j, <icx {I, k). We set 1 < Ui. 

• Let {i,j),{k,l) G Nq be distinct. We write -< {k,l) if i > k, or 
i = k and j < I. The order Cn < ■ ■ ■ < ci, induces a ranking on C, 
using the monomial order -<: < c^^i -< {k,l). 

We call r the ranking on C U U that eliminates U with respect to C, that 
is c?^x < d''*u, for all x G C, m G f/ and /c, k* G Nq. The previous are all 
classical concepts in differential algebra and references for them are [19] and 



Let V := {/i, . . . , /n} be a system of linear differential polynomials in 
D{[/}. We assume that: 

(PI) The order of /j is Oj > 0, z = 1, . . . , ra. So that no fi belongs to D. 
{V2) V contains n distinct polynomials. 

(P3) "P is a nonhomogeneous system. There exist G D and hi homogeneous 
differential polynomial in such that fiiU) = ai — hiiJJ) and, for 

some i G {1, . . . , n}, ^ 0. 

We denote by D[9] the ring of differential operators with coefficients in D. 
There exist differential operators Cij G D[5] such that 



D{f/} = ©[u,-fc I j = l,...,n-l, A;gNo]. 



[23]. 



n-l 
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We denote by IS"! the number of elements of a set S. We call the indetermi- 
nates U a set of parameters. The number of parameters of V equals 

iy{V) := \{j e {l,...,n-l} I ^0 for some ^ e {l,...,n}}|. (1) 

In addition, we assume: 

(P4) iy{V) =n-l. 

Let [P]d{c/} be the differential ideal generated by P in D{f/} (see [23]). 
Our goal is to define differential resultant formulas to compute elements of 
the elimination ideal 

[P]o{u} n P. 

The assumption iy{V) = n — 1 guarantees [P]b{(7} H D 7^ {0} and allows the 
codimension one possibility, see Section 5 and Example 2.1(1). Nevertheless, 
in this paper we also deal with subsystems V of V such that i^(P') 7^ \V'\ — 1 
and the study of the consequences of the relation between iy{V') and \V'\ is 
central to this work. 

Examples 2.1. 1. Let us consider the system V = {/i, /2, /s} in U2} 

with 

/i = ai + ai,i,oUi + + 01,2,1^2,1 + 01,2,2^2,2, 

/2 = 02 + 02,2,2^2,2 + 02,2,3^2,3, (2) 

/s = O3 + 03,i,iMi,i + 03,2,1^2,1 + 03,2,2M2,2- 

We assume that every coefficient aij^k is nonzero. The differential op- 
erators describing this system are 

,1 — Ol,l,o + di i^id, Ci^2 — Ol,2,lC^ + Oi,2,2f^ , 
>^^2,1 = 0, £2,2 = 02,2,2<9^ + 02,2,3<9^, (3) 

A,l = 03,l,l9, £3,2 = 03,2,1<9 + 03,2,2<9^ 

Given V' = {/i,/2}, ^CP') = 2 = \V'\ and it is easily seen that 
[fly f2]o{ui,u2} n D = {0}. Thus we cannot use V to eliminate Ui 
and its derivatives hut we can eliminate U2 and all its derivatives. 
Namely if ID) = ©{ui} and if ai,a2 are differential indeterminates, 
by [27], Algorithm 2 then there exist nonzero differential operators 
Ci,C2 G I}[d] such that £i(oi) + £2(02) = £i(/i) + £2(72) belongs 

to [P']s{«,}nl^{0}. 
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2. Let us consider a system V = {/i = ci+£i /2 = C2+£2,i(wi), /s = 
C3 + £3_i(mi)} m wi/i D = /C{ci, C2, C3} and each fi of nonzero 

order. Observe that 

u{V) = 1 < |P| - 1 = 2. 

By [27], Algorithm 2, there exist nonzero differential operators £1, £2 G 
D[9] such that 

Ri = £i(ci) + £2(02) e [/i, /2]d{«i} n D 

and nonzero 1^2, "^^a G D[9] snc/i that 

R2 = V2{C2) + I?3(C3) e [/2, /3]d{«i} H ©. 

r/iws [PJdIui} n D has codimension greater than one, it is generated at 
least by two differential polynomials. 

Given a nonzero differential operator C = J^kef^o^kd'' G ID>[9], let us 
denote the support of £ by ©(£) = {k E Nq \ ak 0}, and define 

ldeg(£) := min ©(£), deg(£) := max6(£). 

For J = 1, . . . , n — 1, we define the next positive integers, to construct con- 
venient intervals bounding the supports of the differential operators Cij, 

IjiV) := min{oi - deg(£ij) | 7^ 0, z = 1, . . . , n}, 

:= min{ldeg(£ij) | 7^ 0, z = 1, . . . , n}, ^ ' 

Given j G {l,...,n — 1}, observe that, for all i such that £jj 7^ we have 

< ldeg(£,,) < deg(£,,) < o, - 7,(P). (5) 
Therefore, for £jj 7^ the next set of lattice points contains 6(£jj) , 

:=[7/^),o.-7,(^)]nZ. 
Finally, to explain the construction of Section 3, we will use the integer 

n-l 

1{V):=Y.1,{V). (6) 
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Example 2.2. LetV be as in Example 2.1(1). We have 

01 = 2, 6(A,i) = {0,l}, 6(£i,2) = {l,2}, 

02 = 3, ©(£2,1) =0, ©(£2,2) = {2,3}, 

03 = 2, ©(£3,1) = {1}, ©(£3,2) = {1,2}. 

Thus 

l,iV) = 0, 7i(^) = l, 7i(^) = l, 
72(^) = 1, l2iV) = 0, 72(P) = 1, 

and-f{V) = 2. 

3. Differential resultant formulas 

Let us consider a subset PS of dV := {9'^/i | i = 1, . . . , n. A; G Nq} and a 
set of differential indeterminates U G {U} verifying: 

(psl) PS = {d^fi \ke[0, Li\ n Z, Li G No, z = 1, . . . , n}, 

(ps2) PS C D[W] and \U\ = |PS| - 1. 

LetiV:= E7=lO^■ 

Remark 3.1. Particular cases of sets PS and lA verifying (psl) and (ps2) 
were given in [6] and [26] (see also [27]). 

1. In [6], Li = N - 0, and U = {uj^k | A; G [0, iV] n Z, j = 1, . . . , n - 1}. 

2. In [26], Section 3, Li = N — Oi — ^ , where 7 := X]j=i 7i' 

7j := min{7j('P), min{oi | Cij = 0, z = 1, . . . , n}}, 

and U = {uj^k I A; G [0, — 7j — 7] fl Z, j = 1, . . . , n — 1}. 
Observe that both choices coincide 2/7 = 0. 

The coefficient matrix Ai{PS,U) of the differential polynomials in PS as 
polynomials in 3[U] is a |PS| x |PS| matrix. 

Definition 3.2. Given PS andU verifying (psl) and (ps2), we call 

det(A^(PS,W)) 
a differential resultant formula for V . 
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It can be proved as in [27], Proposition 16(1) that det(A1(PS, W)) G 
[V]o{u} nro. Therefore, if det(Al(PS, W)) ^ then it serves for differential 
ehmination of the variables U from V. See Examples 5.3. 

The differential resultant formulas for V given in [6] and [26] are deter- 
minants of matrices with zero columns in many cases. Let PS^ := {d'^hi I 
d'^fi G PS}, the set containing the homogeneous part of the polynomials in 
PS. The coefficient matrix 

C{PS,U) (7) 

of PS'^, as a set of polynomials in D[W], is a submatrix of Ai{FS,U) of size 
I PS I X (I PS I — 1). We assumed that V is a nonhomogeneous system, thus if 
A^(PS,W) has zero columns, those are columns of £(PS,W). 

Remark 3.3. The differential resultant formula for V given in [26] is called 
the linear complete differential resultant ofV and denoted 9CRes(P). With 
PS andU as m Remark 3.1(2), dCRes{V) = det(Al(PS, W)). Observe that, if 
Iji'P) 7^ for some j G {1, . . . ,n — 1}, then the columns of C(PS,U) indexed 
by Uj, . . . , Uj^^ {j>)-i are zero. If^^iV) > 7j for some j G {1, . . . ,n — 1} , then 
the columns of C{PS,U) indexed by Uj^N-^j(r)-'y+i • • • ? '^j,N-'yj-^ are zero. 

If N — Oi — 7(P) > 0, z = 1, . . . , n, the sets of lattice points Ij := [0, — 
Oi — 'y{V)] n Z are non empty. We define the set of differential polynomials 

ps{V):={d''f\keh,i = l,...,n}, (8) 

containing 

n 

L:=5^(iV-o,-7(P) + l) (9) 

1=1 

differential polynomials, in the set V of L — 1 differential indeterminates 

V := K- I k G [7,(P), N - 7 .(P) - 7(P)] n Z, j = 1, . . . , n - 1}. (10) 

Let us assume that ps(P) = {Pi, . . . , Pl\. For i = 1, . . . , n and k G Ij, 

Pl{i,k) '■= d^fi, 
i-1 

k) := ^(iV - - 7(P) + 1) + iV - - 7(P) + 1 - A; G {1, . . . , L}. 

h=l 
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The matrix Ai{V) := M.{ps{V), V) is an L x L matrix. We assume that the 
/th row of Ai{V), / = 1, . . . , L contains the coefficients of Pi as a polynomial 
in D[V], and that the coefficients are written in decreasing order with respect 
to the orderly ranking on U. 

Thus, if — Oj — 7(P) > 0, z = 1, . . . , n, we can define a linear differential 
resultant formula for V, denoted by dFRes(V), and equal to: 



In general, we cannot guarantee that the columns of Ai iV) are nonzero, 
as the next example shows. 

Example 3.4. Let V = {fi, f2, Is} , with oi = 5, 02 = 1, and 03 = 1. Let 

fi = ai + £i,i(ui),/2 = a2 + A,2(M2), /s = 03 + ^,2(^*2), with&{Ci^i) = {1,5} 
and e(£2,2) = ©(^,2) = {0,1}. Then -f{V) = i^{V) = 1 and N - oi - 
7(P) = 1, iV - 02 - 7(P) = iV - 03 - 7(P) = 5. Therefore M{V) can he 
defined hut columns indexed by 1*1^3 and Ui^4 are zero. 

We give next, sufficient conditions on V for M.(V) to have no zero 
columns. Let 5'„_i be the permutation group of{l,...,?T. — 1}. A linear dif- 
ferential system V is called differentially essential if, there exist i G {1, . . . , n} 

and Ti G Sn-i such that 



Observe that, if V is differentially essential then = n — 1 but the 

converse is false. Differentially essential systems of generic, non necessarily 
linear, differential polynomials were defined in [20], Definition 3.3 and (12) 
is a new characterization of this requirement in the case of linear differential 
polynomials. 

Definition 3.5. A linear differential system V is called super essential if, for 

every i G {1, . . . ,n}, there exists Ti G Sn-i verifying (12). 

The notion in Definition 3.5 is introduced for the ffist time and its im- 
plications will be studied further in Section 4. Simultaneously, in [21] the 
notion of rank essential (non necessarily linear) system was introduced, and 
it is equivalent to super essential in the linear case. 



dFResiV) := det{M(V)). 



(11) 




(12) 
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Given a super essential system P, it will be proved that dFRes(V) can 
be defined and that the matrix Ai{V) has no zero columns. For this pur- 
pose, given i G {1, . . . , ra}, for every t G Sn— 1 we define bijections /i^ : 
{l,...,n}\{z}^{l,...,n-l}by 

f r{n-j), j = l,...,z-l, 

\T{n-j + l), j=t + l,...,n. ^ ' 

In particular, for Tj, i = 1, . . . , n as in Definition 3.5, let 

/ii :=/i;^,i = l,...,n. (14) 

Lemma 3.6. Given a super essential system V, N — Oi — '^(V) > 0, i = 
l,...,n. 

Proof. Given i G {1, . . . , n}, 

je{l,...,n}\{i} 

By Definition 3.5 and (14), 7^ 0, j G {1, . . . ,n}\{i} and, by (5), 

Oj — l^^(j){V) > 0. This proves the result. □ 

By Lemma 3.6, if V is super essential then the differential resultant for- 
mula 5FRes(P) can be defined. Furthermore, we prove next that M.{V) 
has no zero columns although we cannot guarantee that c}FRes('P) 7^ as 
Examples 5.3(3) and 6.7 show. 

Given a linear differential polynomial / G D{?7}, 

n-l 

f = a + YVj{uj), with a G D and Vj G ©[9], 
we define Gj{f) := 6(I?j). 

Remark 3.7. Let f G D{t/} 6e linear. If k e Gjif) butk + 1^ Gj{f) then 
k + le &j{df). 

Theorem 3.8. Given a super essential systemV (as in Section 2), the matrix 
Ai{V) has no zero columns. 
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Proof. Equivalently, we will prove that ps('P) is a system of L polynomials 
in the L — 1 (algebraic) indeterminates of the set V in (10). We will prove 
that for every j G {1, . . . , n — 1} 

U;ep.(P)©,(/) = [7/^),iV-7,(^) -7(^)] nz. 

Let us denote 'jjiV), IjiV) and Iji'P): j = 1, ■ ■ ■ ,n — 1 simply by jj, 7^ and 
7_^. respectively, in this proof. 

Given j G {1, . . . , n — 1}, the set {deg(£jj) | Cij 7^ 0, i = 1, . . . , n} is not 
empty because i^(V) = n — 1. Let J(j) G {1, . . . , ra} be such that 

deg(£/(j)j) = min{deg(£ij) | dj 7^ 0, z = 1, . . . , n} (15) 

and define d := deg(£7(j)j). Let /(j) G {1, . . . ,n} be such that 

%=^(i)-deg(%),,)- 

Observe that it may happen that = but not necessarily. We can 
write 

[7 ., iV-7j-7] = [7^., d-lp[d, Ar-o,(,)-7+rf]U[Ar-o/(,)-7+rf+l, iV-%-7]. 

If 7^. = d then the first interval is empty. If /(j) = /(j) then d = oj(j) — 7^- 
and the third interval is empty. 

1. For every k G [d, A^ — 07(j) — 7 + 6?] flZ, since d = deg(£7(j)j), by Remark 
3.7 we have k G &j{d''-^fi^j)). 

2. If J(j) ^ 7(j). Given k e [N - o/q) - 7 + d + 1, iV - 7^- - 7]. Let 
(i := deg(£j(-j.) ,^.) = oj^^) — 7^ and observe that 

^ = ^(i) ~ + 7^. < - o/{i) - 7 + 7^. < - o/{i) - 7 + < ^ 
since 7^. < rf. The previous shows that 

k — d<N — ^j— J — d = N — oj(j-) — 7. 

By Remark 3.7, A; G &j{d''-^fj^-X 
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3. li d 'J.. Observe that [7^., d — l]n(U"^;^(3j(/i)) 7^ because it contains 
7.. For k e [l.,d- 1] n (U"=i6j(/i)), there exists ik E {l,...,n} 
^ich that k e~(5j{fi,). If [7.,ci - l]\(Uti©i(/^)) ^ 0> given k e 
[7^.,rf-l]\(ur=i6,(/.)),let 

fc' = max[7^.,A;-l]n(ur=i6,(/.)). 

Observe that k — k' < d — 7 . and there exists i^/ G {1, . . . , n} such that 
k' e 6,(4,) but k' + l^ 6,(4,). If ^fe' ^ /(j) then 

k-k' < oii^j) -jj<N- Oi^, - 7. (16) 

Otherwise, ik' = and k' e &j{fi(^j^). Since P is super essential, 
the bijection given by (14) is defined and for /' = I^J^j^ij) E 
{1, . . . ,n}\{I{j)}, Cjij 7^ 0. In particular, deg(£//j) < o// — 7^- and 

k - k' < d - < Op - < N - oi(^j) - 7. (17) 
By (16), (17) and Remark 3.7 it is proved that k G &j{d^~^' fi^,). 

□ 

Example 3.9. 1. Let us have a new look at the system of Example 3.4, 

fi = ai+ £i,i(ni) + , 

/2 = a2+ + £2,2(^x2), (18) 

/s = 03+ + £3,2(^2) • 

It is differentially essential, namely for i = 3 the permutation = 
(2 1) e 5*2 verifies (12), £i,r3{2) = £1,1 7^ and £2,r3(i) = £2,2 7^ 0. 
This system is not super essential, for i = 1 we cannot find ti G 5*2 
verifying (12). The subsystem V' = {/2, fs} in D{m2} is super essential 
and 9FRes(P') is the determinant of a Ax 4 matrix while dFR,es{V) is 
the determinant of a matrix 14 x 14. 
2. The system in Example 2.1(1) is super essential, let us construct AiiV) . 
Using the information in Example 2.2, L = ^^=i(-/V — Oj — 7(P) + 1) = 
11 and 

PS{V) = {9^1, 5/1, /l,5V2, 9/2, /2, 93/3, 9/3, /s}, 

V = {^2,5, U2,4., «1,4, U2,3, «1,3, «2,2, «1,2, «2,1, Ul,l, 
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whose elements are arranged in the order indexing the L rows and first 
L — 1 columns of Ai (V) respectively. We show next the columns of the 

matrix AiiV), we denote d^aij^k by a-flk ^'^^ '^''^i ^f^> ^ ^ ^''^^ 
space limitations. Observe that the first L — 1 columns are the columns 
of the matrix C{ps{V),V) in (7). 

U2,5 "2,4 Mi,4 M2,3 ^1,3 



— ai,2,2 ^1,2,1 + 3a]^ 2 2 '^1,1,1 3a]^2i + 3ai2 2 c^i.i.o + Sa^ 

d^fi ai,2,2 ai,2,i + Sa^^^^ 

<9/i ^ ai,2,2 ' 

/i ^ ' 

9^/2—)- 02,2,3 '^2,2,2 + 2 3 2 2 "I" '^2 2 3 

d^f2 ^ 02,2,3 02,2,2 + 02!2,3 

/2 ^ 02,2,3 ' ' 

d^fs 0,3,2,2 03,2,1 + 3oy 2 '^3,1.1 303,2,1 + 3O3 2,2 3O3 J -[^ 

d'^fs 03,2,2 03,2,1 + 2o^^^ 2 03,1,1 

5/3-^0 03,2,2 

/3 ^ 

M2,2 Ml,2 M2,l Ml,l Ml 1 

(2) , ^(3) 0^(1) , oj2) (3) ^(2) , ^(3) ^(3) ^(3) 



,1 



'^1,1,0 



'-'"1,2,1 ' "1,2,2 '-'"'1,1,0 ' '-'"1,1,1 "1,2,1 '-'"1,1,0 ' 1,1,1 
^"1,2,1 + "1,2,2 "1,1,0 "T ^"1,1 1 "1 2,1 ^'^1,1 "T '^1,1 1 "-l 1,0 '^1 

, (1) (i) ' , (1) ' (i) (1) 

Ol,2,l + Oi 2,2 "^1,1,1 '^1,2,1 '^1,1,0 + '^1,1,1 '^1,1,0 Oi 

Ol,2,2 '^1,2,1 

«S,2 

02,2 

02,2,2 

3a(2) , .(3) 3^(2) (3) 

'-'"3,2,1 ^ "3,2,2 '-'"3,1,1 "3,2,1 

^"3,2,1 ^ "3,2,2 ^"3,1,1 "3,2,1 
, (1) (1) 

03,2,1 -r 0^ 2,2 '^3,1,1 %,2,1 

03,2,2 C^3,2,l 



Ol,l,l 


Ol,l,0 


Ol 
















a? 








02 


3,1,1 





af> 


"3,1,1 





af 


3,1,1 





4" 


03,1,1 





O3 
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4. Irredundant systems of linear differential polynomials 



A key fact to eliminate the differential variables in U from the system V 
is that not all the polynomials in V have to be involved in the computation. 
Namely only the polynomials in a super essential subsystem V* of V are 
needed to achieve the elimination. In this section, it is proved that every 
system V contains a super essential subsystem V* and a new characterization 
of differentially essential systems is given, namely they are the system having 
a unique super essential subsystem. 

The linear differential system V is an overdetermined system, in the dif- 
ferential variables U. Recall that we assumed i^('P) = n — 1 = \V\ — 1. It 
is proved in this section that, the super essential condition on V is equiva- 
lent with every proper subsystem V' of V not being overdetermined, in the 
differential variables U. A name for this idea seems to be lacking in the 
literature. 

Definition 4.1. A system of linear differential polynomials V is called irre- 
dundant (for differential elimination purposes), if every proper subsystem V' 
ofV verifies \V'\ < i^(V'). Otherwise, V is called redundant. 

Furthermore, it will be shown in this section that every linear differential 
system V (even if it is not differentially essential) contains a super essential 
subsystem V*. Let Vi:=P\{fi}. 

Proposition 4.2. IfV is super essential then V is irredundant. 

Proof. For every proper subset V = {fhi, ■ ■ ■ , fhm} of V, there exists i G 
{!,..., 77,} such that V C Vi. Therefore hi,...,hm G {1, . . . , n}\{z} and 
given fii as in (14), 



Let Xj j, i — 1, . . . , n, J — 1, . . . , n — 1 be algebraic indeterminates over Q, 
the field of rational numbers. Let X{V) = [Xij) be the n x [n — 1) matrix, 
such that 



We denote by Xi{V), i = l,...,n, the submatrix of X{V) obtained by 
removing its ith row. Thus X{V) is an n x (ri — 1) matrix with entries 
in the field K := Q{Xij \ Xij 7^ 0). 



^ht,Mht) 7^ 0, t = 1, . . . ,m. 
Since fii is a bijection, ^'{V') >m= \V'\. 



□ 




= 0. 



(19) 
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Lemma 4.3. Given i G {1, . . . , n}, det{Xi(V)) 7^ and only if there exists 
Ti E Sn^i verifying (12). 

Proof. Given r G Sn-i, let us consider the bijection /i^ := /i^ as in (13). We 
can write 

det(X,(P))= J2 n ^^-./^^O)- (20) 

rGS„_i je{l,...,n}\{i} 

The entries of XiiV) are either algebraic indeterminates or zero. Thus 
det(Xj('P)) = if and only if every summand of (20) is zero, it contains 
a zero entry. That is, for every r G 5'„_i, there exists j G {1, . . . ,n}\{i} 
such that Xj^^^i^j) = 0, thus Cj^^^^j) = 0. This proves that, det{Xi(V)) = if 
and only if there is no r G Sn-i verifying (12). □ 

Remark 4.4. From Lemma 4-3 we can conclude that: 

1. V is differentially essential rank(X('P)) = n — 1. 

2. V is super essential det{Xi{V)) ^ 0, i = 1, . . . ,n. 

Given the set P := {pi, . . . ,pn} of algebraic polynomials in K[C, t/], K = 
Q{Xij I Xi,, ^ 0), with 

n-l 

Pi Cj + ^ ^ij^j^ i = 1,. . .,n, 
i=i 

a coefficient matrix M(P) of P is an x (2?7. — 1) matrix and it can be obtained 
by concatenating X[V) with the identity matrix of size n. 



M(P) 



1 ••■ 

X{V) 

■■■ 1 



(21) 



The reduced echelon form of M(P) is the coefficient matrix of the reduced 
Grobner basis B = {cq, ei, . . . , e„_i} of the algebraic ideal (P) generated by 
Pi, . . . ,Pn in K[C, U], with respect to lex monomial order with ui > ■ ■ ■ > 
Un-i > ci > ■ ■ ■ > c„, ([8], p. 95, Exercise 10). We assume that cq < ei < 
■ ■ ■ < en-i- 

Observe that the elements of B are linear homogeneous polynomials in 

K[C, f/] and at least 

eo G So := -BnK[C]. (22) 
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Given a linear homogeneous polynomial e G K[C], e = Yl^=iXhChi Xh ^ IK, 
let X(e) := {/i G {1, . . . , n} | X/i 7^ 0}. Let us consider the system 

V* ■.= {h\heX{e^)}. (23) 

Remark 4.5. Lei X := {? G {1, . . . , n} | det(Xi(P)) = 0}. By Remark 44 
the following statements hold 

V is differentially essential -^X ^ {1, . . . , n} and , , 

V is super essential <(=^ X = 0. 

Furthermore, if V is differentially essential then, by Remark 4-4! = {^o} 
and by (24), up to a nonzero constant, 

eo= Yl det(X,(P))Q, with X(eo) = {1, • • • , n}\I, (25) 

iGX(eo) 

the determinant of the matrix obtained by concatenating X{V) with the col- 
umn vector containing Ci, . . . , c„. 

Lemma 4.6. IfV is super essential then V = V* , otherwise V* ^ V. 

Proof. By Remark 4.5, if V is super essential X(eo) = that is 

V* = V. Otherwise, X 7^ and we have two possibilities: if X 7^ {1, . . . , n} 
then, by (25), X(eo) ^ {1, . . . , n}; if X = {1, . . . , n} then, by Remark 4.4(1) 
and (24), rank(X(P)) < n — 1 and Ci G Bq with Cq < Ci, which implies 
X(eo)£{l,...,n}. □ 

We will prove next that V* is a super essential subsystem of V. 

Lemma 4.7. 1. For every V ^ V* , rank(X(P')) = 1^1- 
2. rank(X(P*)) = \V*\ - 1. 

Proof. 1. Given a proper subsystem V' of "P, the matrix X(V') has size 
|P'| X (n — 1). Thus rank(X(P')) < The coefficient matrix in 

K[C', U] of r := {p;, I fh G P'}, with C" = {Ch I A G V'}, is 



M(r) 



1 ■■■ 

X{V') 

■■■ 1 
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If rank(X(P')) < \V'\ then there exists e G (F) n K[C"], the vector 
whose coefficients are in the last row of the reduced echelon form of 
M(r). Therefore, T" ^ V* together with (23) imply 

X(e) C {/i e {1, . . . , n} I /, G P'} g X(eo), (26) 

and e G (P') n K[C] C (P) n = (i3o). This contradicts that 

Bq is a Grobner basis of (P) fl K[C], since rem(e,So), the remainder 
of the division of e by So, equals by (26) rem(e, eo) 7^ 0. Therefore 
rank(X(P')) = [P'V 
2. Let m = \V*\ and P* := {ph \ h e X(eo)}. By 1, rank(X(P')) = m - 1 
for every V ^ with \V'\ = m - 1. Thus rank(X(-p*)) > m - 1 
because X{V') is a submatrix ofX{V*). On the other hand Cq G (P*) fl 
]K[C] implies rank(X(P*)) < m, otherwise the reduced echelon form of 
M(P*) provides no vector in ]K[C]. Therefore rank(X(7^*)) = m — 1. 

□ 

Given a proper subsystem V = {gi := fij^,...,gm '■= of V and 
J = {ji, . . . ,jm~i} C {1, ... ,n — 1}, let F-'('P') be the m x (m — 1) matrix 

y^(p') := (F^ ,), , := X,,^,,^, /i = 1, . . . , m, A: = 1, . . . , m - 1. (27) 

Denote by Yf^iV') the submatrix of Y'^{V') obtained by removing the hth 
row, h = 1, . . . ,m. If P' is super essential then there exists J C {!,..., n—1}, 
\J\ = m — 1 such that: 

9h = flih + XI ^ihji^j)^ /i = 1, . . . , m, (28) 

and 

det(y/(P'))^0, /i = l,...,m. (29) 
That is, i^{V') = \V'\ - 1 and Remark 4.4(2) is verified. 

Theorem 4.8. IfV is not super essential then, the system V* given by (23) 
is a proper super essential subsystem ofV, with v{V*) = \V*\ — 1. 

Proof. We can write V* = {gi := /j^, . . . ,gm '■= fi„,}- By Lemma 4.7, there 
exists J = {ji, . . . ,jm-i} £ ■ ■ ■ ,n — 1}, such that 

det(y^(P*)) 0. (30) 
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Let us denote Y-^{V*) simply by Y(V*) and Y/(P*) by Yh{V*), h=l,...,m, 
in the remaining parts of the proof. Observe that Y{V*) is a submatrix of 
X{V*). We will prove that, the only nonzero entries of X{V*) are the ones 
in the submatrix Y{V*), that is (28) is verified or equivalently 



m—l 



(31) 



k=l 



and also 

det(n(P*)) ^0, /i = l,...,m. (32) 

For this purpose, we will prove the following claims. For / G {1, . . . , m}, 
if det{Yi{V*)) ^ then 



m— 1 



(33) 



k=l 



there exists a bijection r/; : {1, . . . , "m-jyi/} — > {1, . . . , m — 1} such that 

det{Yt{V*)) ^0,\/te Ti := {t G {1, . . . , m}\{/} | ^ 0}. (34) 

1. Proof of (33). Otherwise, there exists j G {1, . . . ,n — 1}\J such that 
Xi^j 7^ 0. This means that the matrix 



X.. 



H,3 



Y{V* 



Xi 



is nonsingular, which contradicts rank(X('P*)) = m — 1, see Lemma 
4.7. 

2. Proof of (34). Since det(y'z(P*)) ^ 0, by Lemma 4.3, there exists 
Ti G ^m-i and a bijection 



r/z:{l,...,m}\{/}^{l,...,m-l} 

ri(m - /i), 

Tiim-h + l), h = l + l 



h = 1,...,/-1, 



such that 



YhMh) 7^0, he {l,...,m}\{/}. 



(35) 
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Given t ^ Ti and the permutation p(/, t) : {1, . . . , m} — > {1, . . . , m}, 
such that 

{t, h = l, 
I, h = t, 
h, he{l,...,m}\{tj}, 

we define the bijection 

?7t : {1, . . . , m}\{t} — > {l,...,m-l},rit = riio p{l, t). (36) 

Thus, by (35) and the definition of Ti, lh,r?t(/i) 7^ 0, /i G {1, . . . , m}\{t}, 
which proves that det{Yt{V*)) ^ 0. 

We are ready to prove (31) and (32). By (30) and (34), aheady (32) holds 
for h e Tm U {m} and, by (33), (31) holds for h G U {m}. We follow the 
next loop to prove (31) and (32) for G {1, . . . , m — l}\Tm. 

1. Set T := Tm and V := {gh \ h \J {m}}. 

2. If T = {1, . . . ,m - 1} then V* = V, which proves (31) and (32), by 
(30), (34) and (33). 

3. liT ^ {1, . . . , m— 1} then, there exists / G T such that T/\(TU{m}) 7^ 
(see below). Set T := {T U Ti)\{m}, V := {gn \ h e T U {m}} and 
observe that by (34), (32) holds for /i G T U {m} and by (33), (31) 
holds for G T U {m}. Go to step 2. 

We prove next that the loop finishes because each time we go to step 3 at 
least one new element is added to T. More precisely, we prove that (in the 
situation of step 3) there exists / G T such that Ti\(T U {m}) 7^ 0. We 
assume that at some iteration Ti C TU {m}, V/ G T to reach a contradiction. 
Given Z G T, if / G T^, by (36), rji = r]m° pijn, I) and 

{Tii{t) \teTi](z {^,,{t) 1 1 G T}, 

else there exist li,...,lp G {l,...,m — 1} such that / G Ti^, Ik G Ti^_^, 
k = 2, . . . ,p and h G T^, by (36) 

{Vi{t)\teTi}C{r^i^{t)\teT\{l,}} = 

W I t e T\{/,_i}} = ■ ■ ■ = {%(t) I t G T\{/i}} C {r/^(t) I t G T}. 

By definition of V and T/, we have proved that i^iV') < \T\ and thus 
rank(X(P')) < i^iV) < \T\, contradicting Lemma 4.7 since V ^ V* and 
|P'| = |T| + 1. □ 
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In particular, Theorem 4.8 shows that if V is not super essential then V 
is redundant, which together with Proposition 4.2 proves the next result. 

Corollary 4.9. A linear differential system V is irredundant if and only if 
it is super essential. 

The next result shows that if V is differentially essential then V* is in 
fact the only super essential subsystem of V. This new characterization of 
differentially essential systems (in the linear case) has now a flavor similar to 
the essential condition in the algebraic case, see [28], Section 1. 

Theorem 4.10. V is differentially essential if and only if V has a unique 
super essential subsystem. 

Proof. 1. If "P is differentially essential, by Remark 4.4 and (22), (P) fl 
K[C] = (eo). By Theorem 4.8, V* is super essential. Let us assume 
that there exists a super essential subsystem V' = {/tj, ■ ■ ■ , fts} of P 
different from V*. This means that {ti, . . . , tg} 7^ X(eo) = {ii, ■ ■ ■ , im}- 
Let r ■.= {p, \ fe V'}, by (29) and (25), (P') n K[C] = (e) with 

s 

e =^det(Y(^(P'))cip and every det{Yi^{V')) ^ 0, 
1=1 

for K C {1, . . . ,n - 1}, \K\ = s - 1. This contradicts that e G (P) fl 
K[C] = (eo) because X(e) = {ti, . . . ,ts} 2r(eo). 
2. Conversely, if V is not differentially essential then, by Remark 4.4, 
rank(X(P)) < n — 1. This implies that the initial variable t G 
{!,..., ra} of Co w.r.t. the order ci > C2 > ■ ■ ■ > c„ verifies l > 2 
because eo is obtained from (21). Let p E Sn he the permutation of 
1 and L. By the same reasoning, if we compute the reduced Grobner 
basis B' = {cq, . . . , e'^_i} of (P) w.r.t. lex monomial order, with 

Ui> ■■ ■ > Un-i > C, > Cp(2) > ■ ■ ■ > Cp(„) 

and Cq < ■ ■ ■ < e^„]^, then the initial variable of Cq is not c^. Thus 
X(eo) 7^ X(eo) and by Theorem 4.8, {/i | i G X(eo)} is also a super 
essential subsystem of V, different from V*. 

□ 
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Examples 4.11. 1. Given the system V = {/i = +£i_2(m2)5 /2 = 

£2,1(^1), /s = '£3^2(^2)} the matrix X{V) defined by (19) equals 



X{V) 



X2,l 
Xgo 



By Remark V is supper essential and, by Lemma 4-6, V* = V. 
2. Let V be a system such that 



X{V) 



By Remark 4-4; is differentially essential but it is not super essential. 
The reduced echelon form of the matrix M(P) in (21) is 





Xl,2 










3^2,3 













X4,2 






E 



Xl,l Xi,2 

X3,2 









3^2,3 





1 



-a;2,i/a;i,i 1 







1 

a;2,ia;i,2/a;i, 1X3,2 

-a;4,2/x3,2 1 



The columns of E are indexed by Ui > U2 > U3 > Ci > C2 > > C4 
and its last row gives the coefficients of cq, see (22). Thus X(eo) = 
{3,4} and V* = {/s, f^} is the only super essential subsystem ofV, by 
Theorems 4-8 and 4-10. 
Let V be a system such that 



X{V) 



Xl,l Xi^2 Xi-^ 

X2,2 

0:3,2 

X42 



By Remark 4-4! "^s not differentially essential and thus it is not super 
essential either. The reduced echelon form of the matrix M(P) in (21) 
is 



1 

-X3,2/X2,2 1 

-a;4,2/x3,2 



E 



xi,i 






a;i,2 
a;2,2 







Xl,3 







1 











1 
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The columns of E are indexed by ui > U2 > > Ci > C2 > C3 > C4 and 
its last two rows give the coefficients of Cq < Ci such that Bq = {cq, ei}, 
see (22). Thus {f3, f^} is a super essential subsystem ofV but in this 
case 1/2, /s} and {/2,/4} are also super essential subsystems ofV. 

5. Differential elimination for systems of linear DPPEs 

In this section, we set D = /C{C} and consider a system of linear differ- 
ential polynomials in D{t/} = /C{C, t/}, 

V = {F, := a - H,{U),i = 1, . . . , n}, (37) 

with -Hi{U) = EJ=i^ A,i(^^i), e IC[d]. Observe that V verifies (P2) 
and (V3) (as in Section 2), let us assume (VI) and (Vi). Let [Plujc/} be 
the differential ideal generated by V in D{f/}. By [13], Lemmas 3.1 and 3.2, 
["^Idic/} is a differential prime ideal whose elimination ideal in D equals 

1D(P) := [V^u} n D = {/ G D I f{H,{U), H^{U)) = 0}. 

It is called in [13] the implicit ideal of the system of linear differential poly- 
nomial parametric equations (linear DPPEs) 

r ci = H,{u). 

\cn = H4U). 

Let PS C dV and U C {U} be sets verifying (psl) and (ps2) (as in Section 3 
but with V as in (37)). The set PS belongs to the polynomial ring /C[Cps,W], 
with 

Cps := {Ci,k I A; G [0, Li] n Z, z = 1, . . . , n}. 

Let (PS) be the algebraic ideal generated by PS in /C[Cps,W]. 

The implicitization of linear DPPEs by differential resultant formulas was 
studied in [27] and [26]. The results in [27] and [26] were written for specific 
choices of PS and U, as described in Remark 3.1. In this section, some of the 
results in [26] are presented for general PS and U, to be used in Section 6, 
namely Theorem 5.2. In addition, the perturbation methods in [26], Section 
6 are extended to be used with formula (9FRes('P). We also emphasize on 
the relation between the implicit ideal of V and the implicit ideals of its 
subsystems. 
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Let V be a subsystem of V. If \V'\ = m then V' = {F^^, . . . , Fh,.^} and 
the imphcit ideal of V equals 

ID(P') = {/ e IC{C'} I f{H,,{U), . . . , H,JU)) = 0}, (38) 

where C = {q | G P'}. Let D' := /C{C"}. If < then it may 

happen that ID(P') = {0}, see Example 2.1(1). 

We use next the notions of characteristic set, generic zero and saturated 
ideal, which are classical in differential algebra and can be found in [23], 
[19], and in the preliminaries of some more recent works as [18] and [15]. If 
\V'\ > i^{V), by [13], Lemma 3.1, ID(P') is a differential prime ideal with 
generic zero {Hh^{U), . . . ,Hh^{U)). Let C be a characteristic set of ID('P') 
(w.r.t. any ranking). The differential dimension of 10(7^') is dim(ID(7'')) = 
m — |C| < m — 1 and coincides with the differential transcendence degree over 
JC of JC{Hh,{U), . . .,HhjU)), see [7], Section 4.2. If V is redundant then, 
there exists V ^ V, with I'iV') < \V'\ and, by the previous observation, 

{0} 7^ ID(P') C ID(P). (39) 

Let U' be the subset of U such that \U'\ = v{V') and V C W{U']. Since V 
is a set of linear differential polynomials, a characteristic set A of [P']©'!!/'}, 
w.r.t. the ranking r, obtained for instance by [18], Algorithm 7.1, is a set of 
linear differential polynomials in ©'{t/'}. If \V'\ > iy{V') then by [13], Theo- 
rem 3.1, ^0 := vAnD' is a characteristic set of ID(P'). By [19], Lemma 2, page 
167 and the fact that the elements in Aq are linear differential polynomials 
in D', 

ID(P') = sat(A) = [A]d', 

where sat(^o) is the saturated ideal of in ID'- If dim(ID(P')) = m — 1 
then ID(P') = [A]oi for a linear differential polynomial v4(c/jj, . . . , Ch^) in D'. 
From the previous discussion we can conclude the following. 

Proposition 5.1. Let V he a proper subsystem of V with I'iV') < \V'\. 
If dim(ID(P)) = n — 1 then ID('P) = [A]n, where A is a nonzero linear 
differential polynomial such that ID{V') = [A]©/. 

Given a nonzero linear differential polynomial B in 10(7^), by [26], Lemma 
4.4 there exist unique J-'i G ]C[d] such that 

n n 

B = Y^Hci) and ^ J-,(i/,([/)) = 0. (40) 

j=l i=l 
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We denote a greatest common left divisor of J^i, . . . , J-'nhj gcld(J-'i, . . . ,J^n)- 
We recall [26], Definition 4.9: 

1. The ID-content of B equals IDcont(i?) := gcld(J^i, . . . , We say 
that B is ID-primitive if IDcont(5) G JC. 

2. There exist £j G IC[d] such that J-^ = IDcont(i?)£j, i = 1, . . . ,n, and 
£i, . . . ,Cn are coprime. An ID-primitive part of B equals 

n 

IDprim(E) := ^^^(q). 

1=1 

If S belongs to (PS) then ord(-B, q) < Lj, i = 1, . . . , n. Given a nonzero 
linear differential polynomial B in (PS), we define the co-order with respect to 
PS of B to be the highest positive integer cps(-B) such that d'^^^^^^B E (PS). 
Observe that, this definition was given in [26], Definition 4.7, for a choice of 
PS. 

Theorem 5.2. Let V be a system of linear DPPEs as in (37). Let PS C dV 
andU C {U} he sets verifying (psl) and (ps2). //dim(ID(P)) = n — 1 then 
ID('P) = [A]n, where A is a linear differential polynomial verifying: 

1. A is ID -primitive and A G (PS) fl ©. 

2. cps(A) = |PS| - 1 - rank(£(PS,W)). 

Proof. We can adapt the proof of [26], Theorem 5.2. □ 

We can also adapt the proof of [27], Theorem 10 (1)<(=^(3) to show that 

det{M{FS,U)) 7^ ^ rank(£(PS, W)) = |PS| - 1. (41) 

As observed in Section 3, if A = det(A^(PS, W)) 7^ then A is an element 
of the differential elimination ideal ID('P) fl©. The next examples illustrate 
this statement in the case of the formula 9FRes(P). 

Examples 5.3. Let JC = Q{t) and d = -§-^. 

1. Let us consider the system V = {Fi, F2, F3, F4} in 3{ui,U2,U3} with 
D = /C{ci, C2, C3, C4} and 

Fl = ci + 5mi,2 + 3m2 + M3, 

F2 = C2 + Ui+ M3, 

F3 = C3 + Ui^2 + U2 + U3, 

Fi = C4^ + Ui+ M2,l + M3,2- 
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Observe that N = 6 and -f(V) = %{V) = 1. Thus dFRes(V) is the 
determinant of the coefficient matrix of 

ps(p) = {a^Fi, . . . , Fi, a^Fa, . . . , F2, a^Fs, . . . , F3, a^Fi, . . . , F4}. 

Namely 

A =aFRes(P) = £i(ci) + £2(02) + ^(cs) + A(c4) ^ 0, 
A =128c4 + 192c3 + 64503 - Udci + 128^204 - 1289^02 + Ud'^Ci 
- 320(9^ C3 + 64(9^ci + 2569^ cs - 192(9^03 - 64ci - 128c2, 

where 

£1 = 64((5-l)(5 + l)2, 
£2 = -128(5 - -5^-5-1), 
£3 = 192 + 645 - 320(5^ - 1925^ 
£4 = 128 + 12851 

Thus £i,£2,£3 and £4 are coprime and A is an ID-primitive linear 
polynomial in ID(P) fl D. 

2. Let X be a differential indeterminate over /C. Set D = IC{x} and 
specialize Ci = x, C2 = C3 = C4 = in the previous system V to 
obtain = {fi, f2, fs, fi} , which is not a system of DPPEs any 
more. Still 9FRes('P^) = GAd'^x + QAd^x — Q4dx — 64x, the special- 
ization ofdFResiV), is an element of the differential elimination ideal 

3. // we replace Fi in V by Ci + U12 + 3^2 + then (9FRes(P) = 0. 

If 9FRes(P) = 0, as in Example 5.3(3), then this formula cannot be used 
to obtain an element of the differential elimination ideal. The perturbation 
methods in [26], Section 6 are next extended to achieve differential elim- 
ination via differential resultant formulas, even in the dY'Res{V) = case. 
Although differential elimination can be always achieved via perturbations of 
9CRes('P), as explained in [26], it is worth to have similar methods available 
for 9FRes('P) (and other possible formulas) since 9FRes('P) is the determi- 
nant of a matrix of smaller size than the matrix used to compute c?CRes(P) 
in many cases ( see Example 5.10). 
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Let p be an algebraic indeterminate over /C, thus d{p) = 0. Denote 
JCp = IC{p) the differential field extension of /C by p. A linear perturbation of 
the system P is a new system 

Ve = {Fr.= F,~pe,{U)\z = l,...,n}, 

where the linear perturbation e = (£:i(f/), . . . , e„(f/)) is a family of linear 
differential polynomials in /C{f/}. The rest of this section is dedicated to 
prove that, if V is super essential then there exists a linear perturbation e 
such that dFResiVe) ^ 0. For this purpose, it is shown how the proof of [26], 
Theorem 6.2 applies to a more general situation than the one in [26]. 

Let us consider cJ := (wi, . . . , aj„) G Nq and /3 := . . . , (in^i) G Nq"^ 
verifying: 

If := YTi=i and /3 := EJ=i then - - /3 > 0, i = 1, . . . , n. 
(/32) If A J 7^ then < deg(£ij) < - 

Taking assumption (7^4) into consideration, for every i there exists j such 
that by (/32) 

< 1] - - /3 < 1] - /3j - /3. 
So it can be easily verified that the next sets satisfy (psl) and (ps2): 

ps/^'^(P) := {d^Fi I G [0, - - /3] n Z, z = 1, . . . , n}, 
U^^^ ■= {ujf, I k G [0,n- (3j - (3]r]Z,j = l,...,n- 1}. 

Thus, we can define the matrix A^(ps'^'^(P), W^'^), as in Section 3. 

Given a nonzero differential operator £ = XlfceNo ^^^'^ ^ ^[d], let := 
afc. Let us define the cJ)-symbol matrix a(j^jj^{V) = (cTij) of V as the 
rax (ra — 1) matrix with 

"^^■^^ 1 0, A, = 0. ^^^^ 

If r2 > 1, we can consider the coefficient matrix Ai{ps^'^ {V)) of 
psI'^iV) := {^'^i/, I A; G [0, fi-Wi-/3-l]nZ, z G {1, . . . , ra}, fi-a;—/?-! > 0}, 

i.e. the submatrix of A^(ps'^'^('P), W^'^) obtained by removing the columns 
indexed by 1 and n-/3j-/35 j = 1, . . . , ra, — 1 and the rows corresponding to 
the coefficients of i = 1, . . . ,n. 
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Remark 5.4. Let 7^ be as in Remark 3.1. In [27] and [26], Section 3, 
Pj = 7j, i = 1, . . . , n — 1, UJi = Oi, i = 1, . . . ,n and, (3 E N^^^ and U E Nq 
verify ([51) and ([52). In fact, 

aCRes(P) = det(A^(ps^'"(P),W'^'^)) and 
dCRes^{Hi, ...Hn) = det(7W(psf'^(P))). 

Using the matrix cr --^{[P] we can adapt the proof of [27], Theorem 10 
{l)^{2) to show that 

det(7W(ps^'^(P),W^'^)) 7^ ^ det(A^(psf'''(P))) ^ 0. (43) 

Let us assume, in addition, that oJ G Nq and /3 E Nq~^ verify 
(/33) Ui - [3n-i > 0, z = 1, . . . , - 1, 

to define the hnear perturbation = (0i(f/), . . . , (pniU)) by 

{Mn-l,a;i-/3„_i, ^ = 1, 

+ M„-i+l, 2 = 2, . . . , n - 1, (44) 
ui, i = n. 

Let (T^_^('P) be the (n — 1) x (n — 1) matrix obtained by removing the nth 

row of (J(p;^){V). 

Lemma 5.5. Let us consider uj E Ng and /3 E Kq"^ verifying (Pl), (P2), 
([i3) and ujn > Wn-i > ■ ■ ■ > i^i- Given the linear perturbation cf) defined by 
(44) it holds that det(aJ^^_^(P)) ^ and det(Al(psJ^(P0))) 7^ 0. 

Proof. The proof of [26], Proposition 6.1 holds under assumptions (/31), (/32), 
(/33) and cj„ > ujn-i > ■ ■ ■ > wi, it does not make use of the precise definition 
of [ij or Ui in [26] (see Remark 5.4). Thus we can adapt the proof of [26], 
Proposition 6.1(2) to conclude that det(A^(psf ^(P^))) ^0. □ 

From (43) and Lemma 5.5 the next result follows. Observe that [26], 
Theorem 6.2 coincides with the next proposition for /3 and uj as in Remark 
5.4. 

Proposition 5.6. Let us considers E Nq and /3 E Ng^"^ verifying ([^1), 
([i2), ([i3) and Un > ojn-i > • • • > i^i- Given a linear system V as in (37), 
there exists a linear perturbation such that 

det(7W(ps'^'^(P^),W^'^))^0. 
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We denote ■jjiV), IjiV) and Iji'P) simply by 7^, ■jj and 'j^ in the remain- 
ing parts of this section. 

Remark 5.7. There exists a system Q = {Gi \ i = 1, . . . ,n} in D{f/} such 
that 

Gi{Ui^y Un-l,^ ) = Fi{Ui, . . . , U„_i). 

—1 —n — 1 

Observe that ord(G'j) < Oi, i = 1, . . . ,n and 7j(Q) may be different from 7^. 
Let uJ = (oi, . . . , On), [5 = (71, . . . , 7n-i) and observe that A^(ps^'^(Q), W^'*^) 
may be different from M.{Q) (as in Section 3) but Ai{ps'^'^{Q),V(^'^) is 
obtained by reorganizing the columns of Ai{V), thus 

dFRes{V) = ±det(A^(ps^'^(Q),W^'^)). 

We may assume w.l.o.g. that o„ > o„_i > ■ ■ ■ > Oi and otherwise rename 
the polynomials in V. Let us assume that V is super essential and let /x := /i„ 
be as in (14), thus 

Oi -7Mi) > 0,2 = l,...,n- 1. 
We define the linear perturbation e = {ei{U), . . . , en{U)) by 

{W^(l),oi -7,^(1), i = ^, 

^A'(0,o.-7,(.) + ^/^(*-i):7^(.„i) ' ^ = 2, . . . , n - 1, (45) 
u^,(n-l),^ , i = n. 

We use the perturbation given by (45) to prove the next result, but there are 
other perturbations that may serve the same purpose. 

Theorem 5.8. Given a super essential system V as in (37), there exists 
a linear perturbation e such that the differential resultant 9FRes(Pe) is a 
nonzero polynomial in /C[p]{C}. 

Proof. Let e be the linear perturbation defined by (45). Let Q = {Gi \ 
i = 1, . . . , 77,} be as in Remark 5.7. U tp = [ipilU), . . . , ipniU)) is the linear 
perturbation defined by 

?^mW,o-7m(>) + i = 2,...,n-l, (46) 

UiM{n-i), i = n 
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then Gf{ui^^^, . . . = F^ui, . . . For (3 = (71, . . . ,7„_i) and 

uj = (oi, . . . , On), by Remark 5.7 

dFResiVe) = ±det(A^(ps^'^(Q^),W^'^)). 

We can assume w.l.o.g. that o„ > o„_i > ■ ■ ■ > oi and that /i = (n — 1, . . . , 1) 
(otherwise rename the elements in U), thus ip = (p, the perturbation given by 
(44). By Lemma 5.5 and (43) det(Al(ps'''^(Q^),^^'^'^)) ^ and the result is 
proved. □ 

By Theorem 5.8, if V is super essential then there exists a linear perturba- 
tion e such that c}FRes(Pe) is a nonzero differential polynomial in /C[p]{C}. 
Let be the lowest degree of (9FRes(P£) in p and A^^ the coefficient of p^^ 
in 5FRes(Pe). It can be proved as in [26], Lemma 6.4 that Ad^ G (ps('P))nD. 
Let As{V) be an ID-primitive part of Ad^. Thus the polynomials Ad^ and 
A^{V) can both be used for differential elimination of the variables U from 
the system V. 

Lemma 5.9. A^iV) is a nonzero linear ID -primitive differential polynomial 
m (ps(P)) n D. 

The next example illustrates these perturbation methods. Computations 
were performed with Maple 15 using our implementation of 9FRes(P). Mak- 
ing this implementation user friendly and publicly available is left as a future 
project. 

Example 5.10. Let V he the system of Example 5.3(3). The perturbed sys- 
tem Ve with e as in (45) and fi = fi^ = (3, 1, 2) is 

Ve = {Fi - pU3^2, F2 - p{Ui + Ug), F3 - ^(^2,1 + Ui), F4 - PU2} ■ 

By Theorem 5.8 we have (9FRes(Pe) ^ 0, in fact 

aFRes(Pe) = P^P{ci, C2, C3, C4), with P = Ad, + pA' e /C[p]{C} 

and Ad, = J^i(ci) + J'2(c2) + J^sics) + J^4{ca), with 

= -24(5 - 1)((5 + 1)2 = C{-6 - 1), 
j-2 = -48(5 - 1)(5 + 1)(52 + 1) = £(-252 _ 2), 
J-3 = 24(5 - 1)(5 + 3)(5 + 1) = C{6 + 3), 
j:; = 48(5-1) (5 + 1) = 2£, 
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and C = 24(6 - + 1). Finally 



AeiV) = -dCi - Ci - 29' C2 - 2C2 + 9C3 + 3C3 + 2C4. (47) 

Observe that dY'Res{Ve) is the determinant of a matrix of size 18 x 18. // 
we use 9CRes(P(^), with as in [26](4), to eliminate the differential vari- 
ables Mi,M2,M3,M4, the polynomial (47) is also obtained but 9CRes(P<^) is the 
determinant of a matrix of size 22 x 22. 



6. Sparse linear differential resultant 

The field Q of rational numbers is a field of constants of the derivation 
d. For z = 1, . . . , n and j = 1, . . . , n — 1, let us consider subsets &ij of No to 
be the supports of differential operators 



= 0, 



whose coefficients are differential indeterminates over Q in the set 

C := Ur=i {c,,,- I k G 

Let Fj, i = l,...,n be a generic sparse linear differential polynomial as 

follows, 

n— 1 rt— 1 

¥i := Q + ^ = Cj + ^ ^ Ci,j,kUj,k- (48) 

In this section, /C = Q(C), a differential field extension of Q with derivation 
d, and D = /C{C}. Consider the system of linear DPPEs in D{f/} 

fp:= {¥i = c^-Mi{U) \i = l,...,n}. 

Let us assume that the order of Fj is Oj > 0, i = 1, . . . , n so that, if Qij ^ 0, 

©„• c /,,(^) = [7/^), - %-(^)] n z. 

We also assume (P4) and observe that (P2) and (P3), in Section 2, are 
verified. 

By [20], Corollary 3.4, the dimension of lD(q}) = ^\^{u} n © is - 1 
if and only if ^ is a differentially essential system. In such case, ID(*P) = 
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sat(i?), the saturation ideal of a unique (up to scaling) irreducible differential 
polynomial i?(ci, . . . , c„) in D = /C{C}. By clearing denominators when 
necessary, we can assume that R G Q{C, C}. By [20], Definition 3.5, R is 
the sparse differential resultant of We will denote it by 9Res(^) and call 
it the sparse linear differential resultant of 

Remark 6.1. Given a differentially essential system by Theorem 5.2, 
ID(^) = [9Res(*P)]D and 9Res(^) is a linear ID-primitive differential poly- 
nomial in ID(^). Observe that 9Res(^) is the implicit equation of the sys- 
tem of linear DPPEs as defined in [27], Definition 2. Furthermore, given 
PS C andU C {f/} verifying (psl) and (ps2), it holds that: 

1. (9Res(<p) belongs to (PS) n © and, 

2. cps(aRes(q3)) = |PS| - 1 - rank(£(PS, W)). 

Proposition 6.2. Let be a differentially essential system. Given PS C 
and U C {U} verifying (psl) and (ps2), the following statements are 
equivalent: 

1. det(7W(PS,W)) ^ 0. 

2. ord(9Res(^), Cj) < Li, i = l,...,n and there exists k G {l,...,n} 
such that ord(9Res(*P), Cfc) = L^. 

Furthermore, if det{M{PS,U)) ^ then det(A4(PS, W)) = a(9Res(q}) for 
some nonzero a G /C. 

Proof. By (41), 1 is equivalent to rank(£(PS, W)) = |PS| — 1. Furthermore, by 
Remark 6.1(2), it is equivalent to cps(9Res(*P)) = and, since 9Res(^) G 
(PS), this is equivalent to 2. Finally, if D = det(A^(PS, W)) ^ then 
D G (PS) n© and Cps(-D) = as well. Since (9Res(^) is ID-primitive, there 
exists a nonzero a G /C such that D = adKes{^). □ 

If ^ is differentially essential then there exists a unique super essential 
subsystem ^* of by Theorem 4.10. If *P is super essential then ^* = 
otherwise, by Theorem 4.8, ^* can be obtained by (23). 

Lemma 6.3. Given i G {1, . . . ,n}, if det{Xi{^)) ^ then, for every subset 
of^i, the differential ideal ID(^') contains no linear differential polyno- 
mial. 
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Proof. Let ^' = {F;,,, . . . , F;,„,}, with h,...,h^ G {1, . . . ,n}\{i}. By 
Lemma 4.3, there exists /ij as in (14), such that Qht,fj.i{ht) 7^ 0, t = 1, . . . , m. 
Given a hnear differential polynomial B G ID(^'), by (40), there exist 
J-,,, . . . , e IC[d] such that Et=i ^/^*(H/.,(t/)) = 0, (5 = E;^! -^/..(c/.J). 
Replacing by zero the coefficients of Qhtj, for t = 1, . . . , m and j 7^ fJ'i{ht), 
this would contradict that . . . are differentially independent. 

This proves that B does not exist. □ 

From the previous lemma, we conclude that to compute 9Res(^) we need 
to use all the elements of 

Corollary 6.4. Let ^ be a differentially essential system with super essential 
subsystem The following statements hold: 

1. (9Res(*P) is a linear ID-primitive differential polynomial in ID(^*). 

2. qj* is the smallest subset of%i such that dRes{^) E ID(q3*). 

3. //9FRes(^*) ^ then dRes{^) = ^dFRes{^*), for a nonzero differ- 
ential polynomial a in K,. 

Proof. By Proposition 5.1, 9Res(^) = 9Res(*P*), which proves 1. Statement 
2 is a consequence of Lemma 6.3 and statement 3 follows from Proposition 
6.2. □ 

Observe that the extraneous factor a does not depend on the variables in 
C = {ci, . . . , Cn}, since a is a nonzero differential polynomial in /C = Q(C). 

Example 6.5. Let us consider the following system ^ in ©{ni, U2} 

Fi = Ci + Ci^i^o^i + 01,2,1^2,1, 
F2 = C2 + 02,1,2^1,2, 

= C3 + C3,i,oUi + C3,2,lM2,l- 

The matrix X{^) is as in Example thus ^ is super essential. The 

formula 9FRes(^) is the determinant of the matrix whose rows can 
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be reorganized to get 



Cl,2,l 





2 9ci,2,i 


Cl,l,0 


<9^Ci,2,l 


2 9ci,i,o 









Cl,l,2 
















9C2 


C3,2,l 





2 9C3,2,1 


C3,l,0 


5^C3,2,1 


2 <9c3,i,o 


5^C3,i,o 










Cl,2,l 





(^Ci,2,l 


Cl,l,0 


f^Ci,i,o 


dci 











Cl,l,2 











C2 








C3,2,l 





<9C3,2,1 


C3,l,0 




9C3 














Cl,2,l 





Cl,l,0 


Cl 














C3,2,l 





C3,l,0 


C3 



iVame/?/ (9FRes(q3) = -ci,i,2<9Res(q3), with 

ord((9Res(^), ci) = 2, ord(«9Res(q3), cs) = anc? ord(«9Res(q3), C3) = 2. 

Let us denote by cr(^) the cj)- symbol matrix of *p as in (42), for 
^= (7i(*?^),---,7„-i(^)) andcJ = (oi,...,o„). 

Lemma 6.6. //rank(o-(q3)) < n - 1 i/ien 9FRes(q3) = 0. 

Proof. Observe that the submatrix of Ai (^) whose columns are indexed by 
^j.A^-7jCP)-7CP)' J = 1) ■ ■ • 5 ''^ ~ 1 has as nonzero rows the rows of o"(^). □ 

Examples of systems ^ with iai\k{a{^)) = n — 1 such that (9FRes(^) = 
have not been found so far, but their non existence has not be proved either. 
In addition observe that, even if c}FRes(*P) 7^ 0, for an specialization "P of ^ 
it may happen that 9FRes(P) = 0, see Example 6.7. 

Example 6.7. Given the differentially essential system of generic differential 
polynomials <:p = {Fi, F2, F3, F4}, 

Fi = Cl + Ci,i,oMl + Ci,i_iMi,i + Ci,3,oM3 + £1,3,1^3,1' 

¥2 = C2 + C2,2,0«2 + C2,2,l«2,l, 

IF3 = C3 + C3,i,oMl + C3,3,oM3, 

F4 = C4 + C4,i,oMl + C4,2,oM2 + £4,3,0^3, 

let us consider the specialization Vof^ 

V = {Ci + Ui + 2ni,i + ^3 + 2^3,1, C2 + ^2 + U2,l, C3 + Ml + U3, C4 + Ml + M2 + ^3}- 

It holds that 9FRes(^) 7^ but c}FRes('P) = 0, even thought V is super 
essential and Ai{V) has no zero columns. We can check, applying [26], 
Algorithm 7.1, that dimlD(P) < 3. 
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By Theorem 5.8, if *P is super essential then there exists a hnear perturba- 
tion e such that 9FRes(*P£) is a nonzero differential polynomial in /C[p]{C}. 
Let Air{^) be as in Section 5. Since Ai^i^) is a nonzero linear ID-primitive 
differential polynomial in (ps(*^)) flD, by Lemma 5.9, the next result follows. 

Proposition 6.8. Given a super essential system there exists a linear 
perturbation e and a nonzero a E K, such that a9Res(*p) = Air{^). 

From Proposition 6.8 we can derive bounds for the order of 9Res(*P) in 
the variables Ci, . . . , c„. Let us consider a differentially essential system ^, of 
generic sparse linear differential polynomials, and the super essential system 
^* of q3. If /* := {z e {1, . . . ,n} I Fi G ^3*} and N* := Y^i^i* Oi then, for 
z = 1, . . . ,n, 

ord(aRes(*P), a) = -1 if z ^ /*, 

ord(9Res(^), a) < N* - o, - 7(^3*) if i e I*. ^ ' 

It was proved in [26] that (9CRes(^<^) 7^ 0, for a linear perturbation 0, see 
Remark 3.3. Thus, given a differentially essential system ^ and reasoning 
as above ord(9Res(*P)) < N — Oi — ^, i = 1, . . . ,n. Observe that (49) is an 
improvement of iV — Oj — 7 whenever ^ is not super essential or if 7(*^3) > 7. 

It cannot be said that (49) are the best bounds in the linear case, since 
an improvement (in some cases) has been just presented in [21], 5.2. But 
at least (49) are the best bounds obtained so far from differential resultant 
formulas. 

7. Concluding remarks 

In this paper, a global approach to differential resultant formulas was 
provided for systems of linear nonhomogeneous ordinary differential poly- 
nomials V. In particular, the formula 9FRes(P) was defined, which is an 
improvement of the existing formulas since, it is the first one given as the 
determinant of a matrix Ai {V) with nonzero columns, for V super essential. 
In addition, every system V was proved to have a super essential subsystem 
V* and therefore the formula dFKes{V*) can be computed in all cases. 

Still 9FRes(P*) may be zero and for this reason, given a system V of linear 
DPPEs, the existence of a linear perturbation e such that 9FRes('P*) 7^ was 
proved. Therefore, to achieve differential elimination of the variables U of 
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the system V the polynomial A^{V*) can always be used, which is a nonzero 
linear ID-primitive differential polynomial in ID('P) fl D, as in Section 5. 

Certainly, there is room for improvement, regarding differential resultant 
formulas in the linear case. For a generic system *p of sparse linear non- 
homogeneous ordinary differential polynomials, 9FRes(*^*) may be zero in 
some cases. Finding a differential resultant formula for ^* that is nonzero in 
all cases would improve further the existing ones. It would be given by the 
determinant of a matrix of smaller size (in some cases) and the candidate to 
be the numerator of a Macaulay style formula for 9Res(*P*). Nevertheless, 
even in such ideal situation perturbation methods will be needed to use such 
formulas for applications, namely to perform differential elimination with 
specializations V of ^, see Example 6.7. 
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